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We examine the quantum dynamics of cold atoms subjected to pairs of closely spaced <5-kicks 
from standing waves of light, and find behaviour quite unlike the well-studied quantum kicked rotor 
(QKR). Recent experiments [Jones et al, Phys. Rev. Lett. 93, 223002 (2004)] identified a regime 
of chaotic, anomalous classical diffusion. We show that the corresponding quantum phase-space has 
a cellular structure, arising from a unitary matrix with oscillating band-width. The corresponding 
eigenstates are exponentially localized, but scale with a fractional power, L ~ ft -0,75 , in contrast 
to the QKR for which L ~ ft -1 . The effect of inter-cell (and intra-cell) transport is investigated by 
studying the spectral fluctuations with both periodic as well as 'open' boundary conditions. 
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I. INTRODUCTION 

A recent experimental study [l| of cesium atoms sub- 
jected to pairs of near <5-kicks using pulsed optical po- 
tentials showed behaviour surprisingly different from the 
well-studied single-kick system, the quantum kicked rotor 
(QKR or 8-KR). The QKR is possibly the most stud- 
ied experimental and theoretical paradigm of classical 
Hamiltonian chaos. The classical dynamics of the sys- 
tem make a gradual transition to chaotic dynamics as a 
function of an effective kick strength K (related to the 
intensity of the optical potential), ft has also been well- 
investigated experimentally, using mainly cesium atoms 
2\. In the large K (chaotic) regime, typical classical 
trajectories are diffusive; to lowest order, the diffusion 
is a random walk in momentum, with diffusion constant 
D w Then, for a corresponding ensemble of atoms, 
say, the average kinetic energy increases linearly with 
time i.e. (p 2 )/2 ~ Dt. Short-ranged correlations to the 
classical diffusion do exist, however, and their effects have 
also been observed experimentally 0] • 

In the corresponding quantum system, the quan- 
tum kicked rotor (QKR), the diffusion is arrested on a 
timescale t* ~ D/ti 2 . In experiments and calculations, 
an initially gaussian momentum distribution N(p,t = 0), 
is seen to evolve into an exponentially localized distribu- 
tion N(p,t > t*) ~ exp — \p\/L, where L ~ j-. This 
quantum suppression of chaotic diffusion is an important 
quantum chaos phenomenon termed Dynamical Localiza- 
tion 0. A formal connection between Dynamical Local- 
ization and Anderson localization, (the exponential local- 
ization of electronic wavefunctions in disordered metals) 
has been made y|. 

Here we consider instead the case where the particles 
are subjected instead to pairs of closely spaced kicks: 
the 2<5-kicked rotor (25-KR). The recent experimental 
study [lj showed that the corresponding classical dif- 
fusion is anomalous, and involves many additional cor- 
rections from weak, but long-ranged (in time) corre- 
lations. The experiment identified periodically spaced 



momentum-trapping regions where the atoms absorb lit- 
tle energy, interspersed by regions where they absorb en- 
ergy even more rapidly than the standard 5-KR system. 
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FIG. 1: Experimental results from Q] showing effects of mo- 
mentum trapping. Each point indicates the energy absorbed 
by a cloud of cesium atoms after dynamical localization. The 
horizontal-axis gives the initial average momenta (p) = Po of 
each cloud of cesium atoms: i.e. each atomic cloud initially 
has a gaussian momentum distribution well-centred on Po. 
The trapping regions correspond to (p 2 ) ~ 0. All the data 
corresponds to effective values of K = 3.2 and h = 1. 



In this work we extend and develop a preliminary study 
of this system @, and reveal unusual localization prop- 
erties and spectral fluctuations not seen in the standard 
QKR. We find a novel dynamical localization regime, 
with generic localization lengths L oc /i -0,75 determined 
by a fractional exponent (~ 0.75). A similar exponent 
was found in P(t) = \(ip(t = 0)\yj(t))\ 2 ~ i -0 ' 75 , the re- 
turn probabilities. To our knowledge, there is no similar 
study of a KAM system so globally characterized by a 
single fractional exponent. A mixed phase-space regime 
(with a mixture of a chaotic regions and stable islands) 
can indeed exhibit fractional return probabilities pj , but 
it will typically have many competing exponents, charac- 
terizing only local regions of phase space. The quantum 
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localization properties are thus not generic, and depend 
sensitively on the detailed phase-space structure. 

For the 25-KR, we present evidence indicating that the 
localization properties are related to the scaling proper- 
ties of phase space in the vicinity of golden-ratio can- 
tori. We argue that the observed behaviour is due to 
the presence of families of golden ratio cantori occupy- 
ing an appreciable fraction of each cell. In addition, the 
'cellular' momentum structure suggested by the exper- 
iments is shown to correspond to a time-evolution op- 
erator with a not previously investigated periodically- 
oscillating banded structure. The corresponding eigen- 
states can be well localized within a single cell defined 
by this oscillating band; or they may extend over several 
cells. 

We characterize the spectral behaviour by two param- 
eters: a filling parameter, R, which quantifies the frac- 
tion of a cell that a typical state occupies, and d, which 
quantifies inter-cell transport. We consider two types 
of distinct (but related) delocalization transitions: (a) a 
— > 1 cell transition as eigenstates fill one cell (R ~ 1 
but d > 2) (b) a 1 — > several cell transition, as states de- 
localize from a single cell to many, occurring for R 3> 1 
and d < 2. 

As the — ► 1 cell transition occurs, all states extend 
into fractal cantori-filled regions bordering each cell; we 
identify a regime with spectral properties (particularly 
spectral variances which, while not identical, may be 
compared with 'critical statistics'. The latter were first 
studied near the critical point of the Anderson Metal- 
Insulator transition (MIT) HH Hind [13. 

So-called 

'critical statistics' are now of much current interest in 
chaotic systems with classical discontinuities (non-KAM 
systems) |T3 . ITM IT3. IT^ : they have been attributed to the 
effect of cantori |lfiul7| : however they are not expected in 
KAM system, due to their non-generic properties. Below, 
we use the term critical statistics in this broader sense, 
rather than specifically the critical point of the MIT. 

Imposing periodic boundary conditions in momentum 
space effectively confines the system's eigenstates to a 
single cell (with toroidal geometry). Such a calculation 
shows a transition from Poisson to Wigner-Dyson statis- 
tics (they are of similar functional form to GOE results, 
so are referred to as such below, although strictly they 
are COE statistics), via a new regime of intermediate 
statistics. We also calculate eigenstates for the open sys- 
tem (non-periodic boundary conditions). This calcula- 
tion yields agreement with the single-cell calculation up 
to the onset of delocalization, but beyond shows rather 
different behaviour. For d ~ 2 the spectral statistics 
show a signature of the onset of delocalization of the 
eigenstates into multiple cells, characterized by a return 
towards Poissonian statistics. 

In summary, our main new results are 1) A study of 
the cellular form of the unitary operator of this system, 
which differs substantially with the behaviour of band 
random matrices (BRM), which approximate the ordi- 
nary QKR. 2) Numerical evidence for behaviour analo- 



gous to critical statistics associated with cantori 3) 
The spectral signatures of delocalization from one cell, to 
several cells. 4) A novel regime of exponential localiza- 
tion determined by the fractional exponent 0.75, which 
coincides closely with the dominant scaling exponent ob- 
tained in the vicinity of golden-ratio cantori. In Section 
II, we review the present state of knowledge of the 2<5-KR. 
In Section III we introduce the time-evolution operator 
and obtain, for the first time, an analytical form for its 
bandwidth. In Section IV we investigate the dynamical 
localization and the fractional exponent v = 0.75. In 
Section V we look at the delocalization within a single 
cell, by solving the problem 'on a torus' in momentum 
space and compare with critical statistics. In Section VI 
we compare the behaviour with a calculation with 'open' 
boundary conditions and the signature for the onset of 
delocalization onto many cells. Finally in Section VII we 
give our conclusions and discussion. 



II. THE 25-KR 

We consider a system with a Hamiltonian correspond- 
ing to a sequence of closely spaced pairs of kicks: 

b 2 

H = y - Kcosx^ [<*(* - nT ) + S (t -nT + e)] , 

n 

where e < T is a short time interval and K is the kick- 
strength. We now use a re-scaled time in units of T. The 
classical map for the 25-KR is then a two-kick map: 

Pj+i = Pj + K sin Xj 
x j+ i = Xj +p j+ ie 
Pj+2 = Pj+i +K sinajj+i 
Xj+2 = Xj+i +p j+2 (l - e). 

(2.1) 

Clearly, the limit e = 1 or corresponds to the Stan- 
dard Map, which describes the classical dynamics of the 
QKR: 

Pi+i = Pi + Ksinxi 

Xi+x = Xi +Pi+i (2.2) 

An experimental realization of this system was ob- 
tained in 0, with cold cesium atoms in pulsed stand- 
ing waves of light. Atoms with momenta (p) = po ~ 
(2m + l)7r/e (relative to the optical lattice) are confined 
in momentum trapping regions and absorb little energy; 
conversely, atoms prepared near p^e = 2rmr experience 
rapid energy growth up to localization. The experimental 
trapping regions are shown in Fig^ The basic mech- 
anism of trapping is fairly intuitive: atoms for which 
Poe = {2m + 1)tt and m = 0, 1, . . . experience an impulse 
Ksvnx followed by another one ~ _ftTsin(a; + 7r) which 
in effect cancels the first. Over time, however, there is 
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a gradual de-phasing of this classical 'anti-resonant' pro- 
cess. A theoretical study of the classical diffusion over 
longer times than a couple of kicks in Q found anomalous 
momentum diffusion for all po, with long-ranged correc- 
tions to the D ~ K 2 /2 uncorrelated diffusion rate, not 
present in the Standard Map. 

Below we present the corresponding quantum be- 
haviour. 



The one-kick matrix for the QKR has a well-studied 
band-structure: since Ji_ m {x) ~ for \l — m\ ^> x, we 
can define a bandwith for U^-°\ namely, 6=4- (this is 
strictly a /iaZ/-bandwidth) which is independent of the 
angular momenta I and m. However, this is not the case 
for the matrix of U e . 

It is shown in the Appendix that, assuming \l — m\ is 
small, we can write: 



III. TIME-EVOLUTION OPERATOR FOR 
25-QKR 

The time evolution operator for this system may be 
written: 



U e — exp —i 



p 2 (T-e) K 



} , ■ exp i— cos x (3-1) 
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FIG. 2: Left (a): Structure of time-evolution matrix for 
a Quantum Kicked Rotor (QKR), in a basis of momentum 
states, showing the constant bandwidth structure typical of a 
Band Random Matrix (BRM). Right (b): Form of U € for our 
system, the 25-QKR, showing the oscillating bandwidth struc- 
ture. Before derealization, eigenstates are confined within a 
single 'momentum cell' of dimension N. 
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(3.4) 



where the phase $ = %[l 2 T + elm + el 2 ] + ir(l - m)/2. 
Hence we infer a momentum dependent bandwidth, 
b(p) = ^jf cospe/2. Fig [21 shows the calculated form of 
both matrices (white denotes matrix elements less than a 
small threshold). While has a constant bandwidth, 
the bandwith for the matrix of U e oscillates with I from 
a maximum value b ma x = ^yf, equivalent to twice the 
bandwidth of EA°) , down to a minimum value 6 m j„ ~ 0. 

In effect, since b m i n ~ 0, U e is partitioned into sub- 
matrices of dimension N = ^n. corresponding precisely 
to the momentum cells of width Ap = Nh observed in 
the experiment. 

For the QKR, the localization properties of the eigen- 
states have been investigated extensively (see e.g. [2£| 
for a review). The eigenstates are exponentially local- 
ized, with momentum probability distributions N(p) ~ 
exp— 2\p\/L, where the localization length L w = 
b 2 H/4:, in the large K, small K limits. 

For the 25-QKR, in the limit of small bandwidth 
2|£ <C N we define a 'local' localization length for the 
eigenstates: 



K 2 



L(p) = b 2 (p)h/4 = — cos 2 pe/2 



(3.5) 



In a basis of plane waves, U e has matrix elements: 

.l 2 h{T-e) 



free rjl — kick 



ut m =ur-u lT 



exp — %- 



A—m 



(3.2) 



where the J n {j-) are integer Bessel functions of the first 
kind. It is easy to see that U 2 ^ klck is invariant if the 
products K e = Ke and h e = he are kept constant; 



while the free propagator U, 



free _ i 



2 simply con- 



tributes a near-random phase. Provided that l 2 Th ^> 2tt, 
the results are quite insensitive to the magnitude of 
(T — e)h. Hence we often find it useful to consider the 
two scaled parameters K e and h t , rather than to vary 
K, e and h independently. 

The result in Eq. (|3.2|) may be compared with the one- 
kick map in Ea l2.2l 



This corresponds quite well with the oscillations seen 
in the experiment in Fig^a); the energy oscillates sinu- 
soidally from a maximum value (p 2 ) ~ 400 ~ 4L 2 ~ AK A 
for K = 3.2, h = 1 , (in re-scaled units) down to a mini- 
mum value (p 2 ) ~ for Pq ~ ir/e. In contrast, Fig^b) 
corresponds to a regime where the eigenstates are tending 
to fill each cell, i.e. L(p) — > Nh. 

In [3 H3 it was shown that the eigenstates of the 
QKR can also be obtained using a matrix of finite 
dimension N. The ratio of the localization length to N, 
was used to characterize the degree of 'filling' of the ma- 
trix E7(°) It was shown that for a ratio R <C 1 the be- 
haviour is Poissonian. With increasing i?, a transition to 
GOE behaviour was observed 

For the 2<5-QKR we can also introduce a 'filling factor' 
R defined by the ratio [l8|: 



K 2 



K 2 



Nh 2 2irh f 



(3.6) 



r/ (0) - exD 



. l 2 Th (K 
~ l 2 l ~ m U 



(3.3) 



While for the 2J-QKR we have a natural choice of N, 
determined by the position of physical boundaries (the 
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trapping regions), the value of N in the QKR case is 
quite arbitrary. The problem is solved 'on a torus' in 
momentum. The momentum periodicity of the matrix 
is adjusted by a choice of a rational value of the kicking 
period Th (see Sec. V below). We can then compare 
matrices for both and U e , with similar N and R. 

The interesting aspect of the 2<5-QKR system is that, 
as we show below, we can vary the coupling between the 
cells independently (with some constraints on allowable 
parameters) from the degree of filling of each individual 
cell. We will show below, that there is a particularly 
interesting regime where, for most eigenstates, L(p) ~ 
NH, but the states are still largely trapped within a single 
cell. We can then 'open' the boundaries of the cells and 
investigate the derealization regime. 

In order to investigate the transport between cells, we 
examine the surprisingly different (relative to the QKR) 
process of dynamical localization in the 2(5-QKR. 



IV. DYNAMICAL LOCALIZATION 
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FIG. 3: Typical momentum distributions, N(p), found after 
dynamical localization, for quantum wavepackets of the 28- 
QKR. Here K = 20, e = 0.0175 and h = 1/8 and 1/30 
respectively. The N(p) of the 25-QKR (slightly smoothed; 
for both eigenstates and wavepackets), show a long-range 
tail of 'staircase' form which on average follows the exponen- 
tial N(p) ~ exp — 2(p — p)/L p where L v /2 — Nh/2d; since 
Nh = 2n/e, the fi-dependence of L p is determined by the 
drop in probability, d, at each cell boundary. 

A set of wavepackets (all initially with N(p,t = 0) ~ 
8(p)) were evolved in time, using the time evolution ma- 
trix U e , for a range of K, e and H. Fig. 10 shows typi- 
cal momentum distributions obtained after a long time 
(beyond the 'break-time' t* for the onset of dynamical 
localization). Fig[3| shows some typical momentum dis- 



tributions N(p). They are modulated by an exponential 
envelope 



N(p) ~ exp - 



2(p - p) 



(4.1) 



but with a regular 'staircase' structure superposed. 
There is a steep drop in probability at each step: 



N+(p) = exp -2d N~(p), 



(4.2) 



where N ± (p) denote the probability after and before the 
step respectively. The localization length is 



_ 2tt 



(4.3) 



The parameter d controls the transport through the cell 
boundaries. It also contains the H dependence of L p . In 
Fig^J the dependence of d on K e and h e — He is shown. 
It may be seen that to very good accuracy: 



d oc 



(fie 



,0.75 



(4.4) 



where f(K e ) is some function of Ke. A rough fit yields 
an estimate, to within 50% or so: 



3.5(fie) 



0.75 



(4.5) 



Hence we obtain the surprising result that for the 25- 
QKR, the localization length has an fi dependence with 
a fractional power, L p ~ K . In comparison, for the 
QKR, i„ ~ fi- 1 . 




Ln(h £ ) 



FIG. 4: Figure shows that Ln(2d) plotted against Ln(he) lies 
close to straight lines of invariant K e = Ke, with constant 
slope w 0.75. Hence d tx (he) ' 75 and L p oc ft -0 ' 75 -in contrast 
to the well-known QKR result L p oc fi -1 . The dashed line 
indicates R ~ 1 the one-cell filling border (i.e. R > 1 below 
the line). The derealization border (d — 2) is the dot-dash 
line (i.e. d < 2 below the line and represents the onset of 
significant coupling between cells). Statistics are presented 
later, in Fig |13l for points corresponding to the dotted line. 

It is interesting to consider the origin of the fi -0 - 75 be- 
haviour. In fact, the first ^1 — 2 steps of the staircase 
can be seen in experiments with optical lattices of pj. 
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An earlier experiment [21| with pairs of broad pulses (as 
opposed to pulses short enough in duration to approx- 
imate 5-kicks) also showed a steep drop in probability 
over a narrow momentum region. This was identified as 
due to the presence of cantori. The finite pulse system 
does not have the momentum periodicity (and hence the 
generic character) of the 28- kicked system, and is effec- 
tively integrable at large p. Nevertheless, the classical 
dynamics is similar to the 2<5-QKR for p ~ 0. The QKR 
shows analogous behaviour: for pulses of finite duration, 
the classical dynamics is similar to the <5-kicked system 
for p small. 

We can examine the classical behaviour in the trap- 
ping regions. Starting from the two- kick map Ea 12. II and 
expanding our initial momenta around the trapping mo- 
menta pj ~ P n + Sp 3 where P n = (2n + l)ir/e, we can 
write 

Pj+2 = Pj + K sin Xj + 
K sin [xj + (2n + l)n + Sp j e + Ke sinxj] . (4.6) 

We can expand the trigonometric expressions, making 
small angle approximations if we assume 8p J e <C 2-7T and 
Ke -C 2n. In the trapping regions we then obtain an 
approximate one-kick map 

Pj+2 — Pj — K 2 -^ sin 2xj — KSp 3 e cos Xj (4.7) 
x j+2 ~ Xj +p j+2 T (4.8) 

The small angle assumption SpP e <C 7r constrains Ap, 
the effective width of the trapping region, 

Ap~l/e. (4.9) 

A study of the Poincare surfaces of section (SOS) may 
be found in |28j |. But in summary, at the centre of the 
trapping region, where Sp° = 0, only the sin 2x kick in 
Ea 14. 71 is significant so the classical SOS shows structure 
very similar to a 'period-doubled' Standard Map with im- 
pulse V'(x) = K e ff sin 2x where K e ff ~ K 2 ^ and K e ff 
is an effective kick strength. However, further out within 
the trapping region, the SOS show that it is the cosx kick 
in Ea 14. 71 which is dominant and the island structure is 
similar to that of a Standard Map with impulse: 

V'(x) ~ KSp ecosx. (4.10) 

Given the importance of the Ke scaling seen in 
FigEl (rather than a K 2 e scaling) we suggest that the 
K e f f cos x form of Eg 14. 71 determines the quantum be- 
haviour. This implies that regions with: 

KSp>e->K 2 ^ (4.11) 

dominate transport. This will be the case if Spj is large; 
from Eo 14.91 we deduce that this implies the criterion 
Ke <C 1. However, if Ke is too small, the phase space 
will be too regular. Hence, values of Ke w 0.1 — 0.3 seem 
indicated. 




FIG. 5: Return probabilities for wavepackets prepared near 
the centre of the trapping regions indicating a Pit) ~ t -0 ' 75 
decay, before localization. The dashed straight lines indicate 
a slope = 0.75. Note that at the lowest value of Ke = 0.25, 
wavepackets prepared at the centre of the trapping region do 
not spread in momentum, (a) h=l/8 (b) h — 1/16 



In any case, the resonance structure in much of the 
trapping regions is always locally similar to the Stan- 
dard Map, but with a varying effective kick strength (and 
phase of the impulse) . For the Standard Map the last in- 
variant phase-space manifolds correspond to p/(2n) ~ <fi 
or ~ 1 — where 4> is the Golden Ratio. The fractal rem- 
nant of this last manifold plays a role in transport in the 
Standard Map for K ~ 1 — 3. Here, this would suggest 
SpKe « 1 — 3 and dp ~ 10. 

The scaling properties of the phase-space around the 
golden ratio cantori were investigated by j2^|. A char- 
acteristic exponent a s w 0.75 was found in directions 
dominated by elliptic fixed points. We term this the sta- 
ble/dominant exponent. On the other hand a u w .66 was 
found in directions dominated by hyperbolic fixed points. 
We term this the unstable/sub-dominant exponent. 

Previous studies of the QKR |23J found a L oc h +a 
scaling, with a ~ a u ~ 0.66, for momenta near the range 
1 — cj) < pT /2ir < <p. In |2]|, a L oc h +<7u scaling was asso- 
ciated to a tunneling type process (termed 'retunneling') 
in |25j. An abrupt change to a regime with L oc h~ a , 
with a w 0.5, was also observed, and was attributed to 
a localization process which dominates when transport 
through the fractal cantoral regions is more 'open'. In 
that work, it was argued that the process has similarities 
with dynamical localization. 

The 2<5-KR results, with a negative sign on the expo- 
nent L ~ ft~ - 75 correspond to a localization regime, and 
are consistent with quantum probability 'sticking' mostly 
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to directions where there are (or were, at lower K) ellip- 
tic fixed points. One might speculate why no previous 
studies uncovered a regime dominated by cr s ; we note 
that in IH 13 the L - h+ - 66 scaling refers to a local 
region; these are regimes where there are many stable 
islands and overall diffusion in each phase-space mani- 
fold is either absent, or rather slow. Here, in contrast, 
part of phase-space is taken by very chaotic, fast diffus- 
ing regions, with trajectories which roam freely over large 
areas of phase-space; the only appreciable localization oc- 
curs at the stabler parts of cantoral remnants. In many 
regimes, there can be contributions from ~ 1/e groups of 
golden-ratio cantori, rather than simply a single region 
near p/2ir ~ <fi like the Standard Map. We note that 
Fig. 13 in |2^|, where the a — 0.5 scaling is obtained, is 
in a difficult numerical regime: barring one data point, a 
value a ~ 0.75 is not implausible. 

We also calculated return probabilities P(t) of the 
quantum wavepackets with time, averaged over 100 ini- 
tial starting conditions close to the centre of the trapping 
region, 

h 

P(t) = (1/100) £<P,(<)>, (4.12) 
h 

where P t (t) = \ (tp(t)\ipi (t = 0))| 2 . The initial condition 
was taken to be an angular momentum eigenstate ipi (t = 
0) = \l). The results were averaged from l\ = n/eh to 
h = h + 100. Fig|5]shows plots of P(t), which show that 
P(t) decays as t~° up to the break-time. For the lowest 
values of Ke, the i -0 - 75 decay is not apparent since the 
wavepacket localizes almost immediately. 

A well-known relation between power-law return prob- 
abilities and the spectral statistics (variances) has been 
investigated for 'critical statistics' in non-KAM billiards 
or the Anderson transition 0, E3, 13 II3- While the 
present systems has important differences (it is a smooth 
KAM system; it has a cellular structure and an oscillat- 
ing band unitary matrix), the pre-eminence of the sin- 
gle fractional exponent motivates an investigation of the 
spectral properties. 

We investigate the statistics as a function of the filling 
factor R and the inter-cell transport parameter d. Calcu- 
lations were done for two types of boundary conditions: 
periodic boundary conditions and 'open' (non-periodic) 
boundary conditions. These are discussed in turn below. 

V. EIGENSTATES: PERIODIC BOUNDARY 
CONDITIONS 

We may keep the dimension of the unitary matrix U to 
a finite value N by making the momentum periodic with 
period NH, following the approach in [lil IT^. |2C| . In 
order to preserve unitarity we must use a resonant value 
of he. We take: 




o is 2 3 

FIG. 6: Typical nearest-neighbour distributions P(s), for pe- 
riodic boundary conditions, for the 25 — KR as a function 
of the filling factor R. The distributions are Poissonian for 
small R, intermediate for R ~ 1 and GOE for sufficiently large 
R. Although the intermediate case never exactly follows the 
semi-Poisson distribution (P(s) = 4se _2s ), a comparison is 
useful, and so is shown in the middle figure 

and then 

Ht = h(T - e) = Mhe (5.2) 

where M is the closest integer to (T — e)/e- We then 
construct a unitary matrix with elements: 

U ln(Tj) = jj exp-iT j {l + 6 x ) 2 /2 

Nl K 

exp i— cos 27r(fc + 9 p ) exp 2ni{k + 9 p )(l — n) 

k=-Ni 

(5.3) 

We may have t\ = he or T2 — mri . Here N\ = (N — 
l)/2. We construct the full two- kick matrix of dimension 
N 

U e (N) = U(tl) U(t 2 ), (5.4) 



(5.1) 



which is then diagonalized to obtain N eigenvalues and 
eigenstates. We may compare the results with the QKR 
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FIG. 7: Comparison of the transition from Poissonian (Q = 0) 
to GOE (Q = 1) NNS statistics for the QKR relative to the 
28- KR as a function of the filling factor R = K 2 /Nh 2 . While 
the QKR makes a rapid transition from Poisson to GOE for 
R ~ 1, the 28 — KR on the other hand, shows a clear change 
in slope where derealization of the quantum eigenstates is 
hindered by the fractal cell borders. Here M = 25; the values 
shown in brackets are approximate values of h = 2irM/N; to 
avoid commensurability of M and N, actual values of N used 
were 629, 1259 and 2513. 



equivalent 



U°(N) = U{TH). 



(5.5) 



Here h = 2Mtt /N ', where M is an integer (non commen- 
surate with N) which determines the momentum width 
of the matrix, i.e. Nh = 2Mtt; unlike the 25-KR, there 
is no underlying physical cellular phase-space structure 
to justify the choice of a particular M, N. 

Note the dependence of the matrix elements on two 
Bloch phases: 9 P is a Bloch phase in the p direction 
while 9 X , the Bloch phase in the x direction, is the quasi- 
momentum. If 9 p ^ 0, n, parity is not a good quantum 
number and we may use all eigenvalues in the statistics, 
regardless of parity. It is customary (e.g. 0,0) to cal- 
culate spectra for several 9 p to improve significance. In 
practice, we found that for the 25 — KR, even for sizable 
K, there are very localized states which are too removed 
from the boundaries and generate pairs of parity-related 
pairs of near-degeneracies. The parity conservation effect 
is more effectively eliminated by, in addition to 9 P ^ 0, it, 
also having 6 X ^ 0, 0.5. The cantori effects we investigate 
are not affected by a non-zero quasi-momentum, and thus 
we average over several quasi-momenta. 




FIG. 8: Typical number variances E2(Z/) as a function of 
R for the 28 — KR. We use the best fit to the slope to 
estimate the level compressibility \. For R small, we have 
Poissonian statistics and ^(L) ~ x L — L. For sufficiently 
large L, E2(L) ~ Ln(L) is far from linear; in this case the 
value x < 0.05 is simply an indication that E2(L) is close to 
the GOE limit. We identify an intermediate regime where 
E 2 (i) ~ xL for 1 « L « iV,(here N = 2513) with x ~ 0.13, 
corresponding to a similar range of R as the regime seen in 

Figia 



A. Effects of cantori in — > 1 cell transition: QKR 
versus 28 — KR 

We have calculated eigenvalues of the unitary matrix 
for a range of values of M, TV and K. Below we compare 
nearest neighbour (P(S)) distributions and spectral vari- 
ances S 2 (i) = {L 2 ) - (L) 2 , of the 25 - KR with those of 
the QKR as a function of the filling factor R = K 2 /Nh 2 . 
We present results for M = 25 (i.e. e = 0.04) and values 
dose to N — 625, 1250 and 2500 (one cannot have N, M 
commensurate) . 

In FigElwe show typical NNS distributions for R ^ 1, 
R ~ 1 and R 3> 1. For both single/double-kicked sys- 
tems, the behaviour evolves from Poissonian to GOE, 
via an intermediate distribution at R ~ 1. While the 
25 — KR does not exactly show the intermediate Semi- 
Poisson form, it can get quite close, and so a comparison 
is helpful. Here we quantify the deviation of P(S) from 
Pp{S) and Pgoe(S), its Poisson and GOE limits respec- 
tively, with a quantity Q |26| : 



l-Q = 



J So (P(S')-P GOE (S')dS') 
J So (P P (S')-P G oE(S'))dS'' 



(5.6) 



Hence Q = indicates a Poisson distribution, while 
Q = 1 signals a GOE distribution. We take So = 0.3. 
In Fig[7] we plot Q as a function of R. We find that 
while the QKR moves rapidly from Poisson to GOE for 
R ~ 1, the 25 — KR curve abruptly changes slope and 
the distribution evolves much more slowly towards GOE. 
We identify this as the regime where delocalization of the 
eigenstates over the whole cell is constrained by the can- 
tori regions bordering the cells. A fit to alternative func- 
tions, such as the Berry- Robnik j^], gives a qualitatively 
similar picture. 

In FigOO we show that the variances are close to Pois- 
sonian ^{L) ~ L for small R, and are close to GOE 
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FIG. 9: Shows slopes (x) of the variances ^(L). Figure 
compares the transition from Poissonian [\ = 1) to GOE 
(X < 0.05) for the QKR with those of the 28 - KR as a 
function of the filling factor R — K 2 /Nfi 2 . Parameters are 
similar to FijQ The QKR makes a rapid transition to GOE 
at R ~ 1 (with the exception of the small minimum seen at 
R ~ 1 and N = 2500 which corresponds to a classical anti- 
resonance K ~ 7r). The 25 — KR on the other hand shows 
an approximate plateau near x — 0.125, the value one might 
expect from the 'critical' statistics relation \ ~ 1/2(1 — D%) 
if D 2 ~ 0.75. 



for large i?. However, for i? ~ 1 there is a regime, with 
nearly linear slope, X^L) — X-^ for L 3> 1 but L <C iV. 
We fitted the slopes of the £2 (L) to the best straight line 
in the range L = 7 — > 37. We note that the graphs are 
not necessarily very linear everywhere: in certain regimes 
there is a pronounced curvature. Nevertheless, the proce- 
dure does give an indication of the average slope. While 
the N = 2500 results may remain linear out to L « 100, 
the N = 625 results saturate at much lower L, hence the 
L = 7 — 37 range is a compromise, good for all three cell 
sizes. 

In Fig|5]we plot the slopes x (the level compressibili- 
ties), for the 2<5— KR. We see that above R ~ 1, there is, 
relative to the QKR, a plateau in the level compressibil- 
ity, where x ~ 0.12. This behaviour is completely absent 
in the QKR: for R ~ 1, the values of x evolve rapidly 
towards the GOE limit. 

We recall the results || for the level compressibilities 
of the Anderson MIT. It is predicted that asymptotically 
the variances have a linear form, for L>1: 




0.001 



100 
KICKS 



1000 



FIG. 10: Return probabilities for the finite dimensional prob- 
lem, for similar parameters to Fig. 5. The P(t) still decay as 
t-°' 75 for t < t*, but P(t) -» 2/N for t -» 00. 



where 



X w 1/2(1 - D 2 /D) < 1 



(5.8) 



and D is the spatial dimension. D2 is a multifractal ex- 
ponent related to the inverse participation ratio. This 
behaviour corresponds Q, in the MIT, to return proba- 
bilities which decay as P(t) ~ t~ D2 and is considered a 
'fingerprint' of critical systems, even in systems without 
disorder such as the non- KAM billiards hj. 

Here we calculate the return probabilities using the 
eigenstates \a) and eigenvalues e iWcv of the N-dimensional 
matrix. We calculate 



I 



-\(l\a)f 



(5.9) 



£ 2 (L)«X£, 



(5.7) 



and average Pi (t) over different starting conditions in the 
trapping regions |Z), exactly as in Eg 14. 121 

We note that in Eg 14. 121 the P(t) were obtained inde- 
pendently of the eigenstates by direct time evolution of 
a wavepacket in the unbounded system. The behaviour 
in Fig^H is very similar to that in Fig|SJ and again a 
decay rate P(t) ~ i~ a75 is apparent. However, P{t) in 
Fig^l] saturates to a larger value than for Fig0 where 
the wavepackets are not restricted to a finite cell. For 
the finite case P(t) — > 2/N as t — ► 00 (the factor of 2 
is attributed to weak localization [2^|). We also observe 
that for R ~ 1, in the regime with x ~ 0.125, wavepack- 
ets started in the trapping regions localize too rapidly 
to demonstrate the power-law decay. Nevertheless the 
behaviour of x is -numerically- consistent with a fractal 
exponent D2 — 0.75. 

We see here that x ~ 0.125 = 1/2(1 — v) would imply 
v « 0.75, close to the exponent we obtained previously. 
This indicates that here v, the dominant exponent near 
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2800 
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<P> 



3400 



FIG. 11: Figure shows the effect of derealization onto multi- 
ple cells: each point indicates the total probability of finding 
a given eigenstate in the cell which contains its average mo- 
mentum (p). Well below the derealization border, (d 3> 2) 
all states are fully contained in their assigned cell (barring a 
few edge-states at the boundary of each cell). But, at the 
'derealization border' (d ~ 2), a proportion of states have 
appreciable probability in neighbouring cells - in fact about 
2% of states have less than half of their probability in their 
assigned cell. Results correspond to e = 0.04,/i = 1/8 and 
K = 3 (upper figure) or K = 7 (lower figure). 



the golden-ratio cantori, may play a role equivalent to the 
multi-fractal exponent in the Anderson MIT. This is 
a significant result as such behaviour has not been seen 
and is not expected in a KAM system [lH ]. 

Unlike the billiard systems, here phase-space is not 
homogeneously filled with cantori. However, underly- 
ing classical trajectories may spend considerable time 
trapped within the fractal trapping regions which bor- 
der the cell. The corresponding typical quantum states 
for R > 1 also sample the cantori region, but localize 
inhomogeneously in the fractal regions. As R increases, 
however, the distribution delocalizes into a more typically 
crgodic regime and the variances evolve towards the GOE 
(£ 2 (L) ~LnL) limit. 



VI. OPEN BOUNDARY CONDITIONS: 
DELOCALIZATION FROM SINGLE TO 
MULTIPLE CELLS 

We diagonalize the matrix of U t , with elements given 
by Ea l3.2l in a basis of a given parity, dimension Ntot = 
10, 000 3> N, for various K,e and fi and obtain all eigen- 



20 
15 ■ 
10 ■ 
5 ■ 



Periodic BC 
open BC 
Poisson 
GOE 



R=0.46 d>2 




R=5 d<2 
K=10 



FIG. 12: Shows the effects of derealization onto multiple cells 
on the spectral variances. Variances obtained with periodic 
boundary conditions (TV ~ 1250) are compared with 'open' 
boundary conditions (7V tot = 10,000 , N « 1250,M = 51). 
For d > 2, the eigenstates are well confined within a single 
cell and there is no significant difference between 'open' and 
'closed' BCs. When d < 2, the results are sensitive to the 
BCs. While the 'closed' (periodic) BCs results tend to GOE, 
the 'open' results revert towards Poissonian statistics. 



values and eigenvectors. For our calculations the cell di- 
mensions are in the range N ~ 600 — 2500 and hence 
each diagonalization spans Ntot/N cells (fewer are kept 
since states in the cells at the edges of the matrix are dis- 
carded) . 

For instance results presented below have e = 0.04 and 
h = 1/8. Hence, each cell contains N = 2n/{he) ~ 1257 
states so each diagonalization obtains up to 7 complete 
cells. Each full spectrum of length Ntot = 10, 000 is 
now split into 8 single cell sub-spectra. We assign the 
i — th eigenvalue to the m— th cell if 



(2m-l)-<pi< (2m+l)-. 



(6.1) 



We calculate P(S) and variances separately for each of 
these sub-spectra of length N = 1257. We then av- 
erage the statistical distributions of 10-40 cells to ob- 
tain smoother distributions. We diagonalize U e for basis 
states between Z = 10, 000 — 20, 000 as easily as states be- 
tween 70, 000 — 80, 000, for example, so as large a number 
of sets of 7 sub-spectra may be obtained as required. For 
the most delocalized spectra (with K = 14)only a sin- 
gle central cell of about 1250 states was sufficiently well 
converged to be used for statistics. 

Fig llOl investigates how much of the momentum prob- 
ability for each eigenstate is contained in the cell it was 
assigned to: for R < 1 and d > 2 the eigenstates are 
essentially fully contained within the cells they are as- 
signed to. However, at the onset of derealization, this 
procedure begins to fail; the expectation value pi of an 
eigenstate may assign it to the m— th cell; however, most 
of its probability may in fact be trapped in neighbouring 
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FIG. 13: Upper figure: shows that the variances (E2 statis- 
tics) for open BCs can pass close to the so-called 'criti- 
cal' regime of near-linear variances, with a slope close to 
X ~ 0.125. Lower figure: compares values of % f° r 'open' 
BCs (a unitary matrix of dimension Ntot = 10, 000 is diago- 
nalized to obtain cells of dimension TV « 1250) with 'closed' 
BCs (a finite matrix of dimension TV ~ 1250, M — 51 is diag- 
onalized, with periodic BCs to retain unitarity). The results 
are in excellent agreement before derealization d > 2. Af- 
ter derealization d < 2, the 'closed' BCs tend to GOE, the 
'open' BCs tend to Poissonian behaviour, providing a clear 
signature of derealization into multiple cells. 



cells. This means that increasingly, the eigenvalues of 
the sub-spectra become uncorrelated and we can expect 
to see a return towards Poissonian statistics. 

The apparent 'failure' of this procedure, in fact pro- 
vides a good marker of derealization from one-to-several 
cells. In contrast, the periodic boundary conditions pre- 
sented in the previous section, could indicate only the 
degree of 'filling' of each single cell. Finally, we compare 
the spectral fluctuations for the two different types of 
boundary conditions (periodic vs open) in Figs ll2lT"ll In 
Fig^Jand Fig^|we show that the spectral variances are 
insensitive to the boundary conditions for d > 2. How- 
ever for d < 2, the variances with periodic boundary con- 
ditions gradually make a transition to GOE behaviour, 
while the variances for the 'open' system begin to return 
back to Poissonian statistics. Fig 1 141 shows that the NNS 
distributions closely follow the same trends. 



VII. CONCLUSIONS AND DISCUSSION 

In conclusion, we have investigated the quantum be- 
haviour of atoms exposed to pairs of 5-kicks and shown 
that the cellular structure arises from a novel oscillatory 
band structure of the corresponding unitary matrix. One 
consequence is a new type of localization-delocalization 
transition not seen in the QKR, where states delocalize 
from a single cell to many, associated with a character- 
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FIG. 14: Behaviour of NNS distributions equivalent to the 
results in Figs ll2l and Fig |13l The results are insensitive to 
BCs before the onset of derealization onto multiple cells d > 
2. After derealization d < 2, the 'closed' BCs yield GOE, 
the 'open' BCs tend to more Poissonian behaviour. 



istic spectral signature (the return to Poissonian statis- 
tics). 

We have also found certain scalings L ~ h~ - 75 and 
P(t) ~ t~ - 75 which we argue result from the bands of 
cantori present on the borders of each cell. We argue also 
that the 0.75 exponent may be identified with the domi- 
nant stability exponent obtained previously near golden- 
ratio cantori. We show that the spectral fluctuations 
(both the NNS and spectral variances) show important 
differences with the QKR in regimes where the dereal- 
ization of eigenstates is hindered by cantori. The numer- 
ics provides evidence for behaviour somewhat analogous 
to 'critical' statistics of non-KAM billiards (in particular 
linear variances, with \ w 0.125). These results are novel 
but important questions remain: while the L ~ /i~ ' 75 
dependence might be explained thus, it remains unclear 
why a similar exponent should also be found in the decay 
of the return probabilities P(t) ~ t -0 75 . The question 
of whether the v — 0.75 exponent is somehow equiva- 
lent to D 2 in the MIT and non-KAM billiards (in other 
words whether v related to any underlying multifractal 
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character of the wavefunctions) remains to be addressed 
in future. 

By implication, the study shows that the behaviour 
of cold atoms in double-pulsed standing waves of light 
is quite different from the single-pulsed systems. Some 
aspects were already identified in the experiments of 1] 
and may have applications in atom optics and atom chips, 
possibly as a mechanism for selecting atoms according to 
their momentum. 
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APPENDIX A: DERIVATION OF THE 
APPROXIMATION fEQlOl) FOR U[ m 

In this Appendix Eg 13.41 will be derived for small e. 
The crucial point in this derivation is that contributions 
to the sum in Eg 13. 21 of terms where 



\k-l\^>K/H \k~m\^>K/h 



(Al) 



are negligible. This implies that for terms that contribute 
appreciably, fc is close to I and to. Therefore the corre- 
sponding operator k = p/h can be written in the form 



k = ko + k\ 



(A2) 



where fco is a number close to I and to. Appreciable 
contributions to the sum are found only for 



| < K/h. 



(A3) 



If he(K/h) 2 <C 1, as is the case for small e, the condition 
hek 2 <C 1 is satisfied, justifying the approximation 

.p 2 e hek 2 hek 2 ~ 

exp— z-^- =exp— i — — w exp — i — — exp — inekoki. 

(A4) 

Substitution in Eq.3.1 with fci = —i-§^i and making use 
of the fact that exp —ihek -^ is the shift operator, one 
finds 



U^ne-iMT-e^+d® (A5) 

i^[cos a;+cos (x-heko)] g-ihekok! ^ 



that reduces to 



U, 



Ira 



,-4[(T-e)/ 2 +efc2+2efc m] 



(A6) 



(l\e 



—i^jr [cos (x—^heko) cos ( ^ heko)] 



HI ) 



The matrix element is calculated with the help of the 
identity 



1 

2^ 



dxe- ll3cos (*-<*) e imx = e lm{a -* /2) J m {f3) (A7) 



resulting in 



u;, 



e '2 



if [(T-e)l 2 +ekl+ek m+ek l] 



(A8) 



^ {l - m) J m -i ^cos(i&fe ) 



Indeed appreciable contributions are found only for | m — 
l\ < 2K/H in agreement with IjAlfl . Since fco m close to 
I and to, within the approximation of this Appendix, it 
can be replaced by one of these. The substitution fco = I 
results in Eo 13.41 
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